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Maximally entangled states are of utmost importance to quantum communication, dense coding, 
and quantum teleportation. With a trapped ion placed inside a high finesse optical cavity, interacting 
with field of an external laser and quantized cavity field, a scheme to generate a maximally entangled 
three qubit GHZ state, is proposed. The dynamics of tripartite entanglement is investigated, using 
negativity as an entanglement measure and linear entropy as a measure of mixedness of a state. 
It is found that (a) the number of modes available to the subsystem determines the maximum 
entanglement of a subsystem, b) at entanglement maxima and minima, linear entropy and negativity 
uniquely determine the nature of state, but the two measures do not induce the same ordering of 
states, and c) for a special choice of system parameters maximally entangled tripartite two mode 
GHZ state is generated. The scheme presented for GHZ state generation is a single step process 
and is reduction free. 

PACS: 03.67.-a, 42.50.-p, 03.67.Dd 



Maximally entangled states are of utmost importance to quantum communication , dense coding , and quantum 
teleportation Q . It is now possible to prepare a cold trapped ion in a given initial state Q with the trap placed inside 
a high finesse optical cavity [f| . We propose an experimental scheme, using ion trap in an optical cavity, to generate a 
maximally entangled three qubit GHZ state. For different sets of initial states of the system, entanglement dynamics 
of ionic internal state, ionic center of mass state, and cavity state is investigated using negativity as an entanglement 
measure and linear entropy as a measure of state purity. 



I. TRIPARTITE SYSTEM AND INTERACTION HAMILTONIAN 



Consider a trapped two-level ion in an optical cavity interacting with an external laser and the quantized cavity 
field. The quantum state of tripartite system (ABC) contains information about the internal state of the ion (sub- 
system A), the vibrational state of ionic center of mass in trap (subsystem B) and the state of optical cavity (C). 
Basis vectors spanning the Hilbert space of composite tripartite system are \i,m,n >, where \i >,i G g, e represents 
ion in ground state (i = g) and excited state (i = e). The vibrational (photonic) number states are denoted by 
\m > (|n >), m(n) — 0, 1, 2, ...oo. The Hamiltonian due to interaction of trapped two-level ion of internal frequency 
luq, with resonant external laser field of frequency u>l = u)q, and with the cavity field tuned to red sideband of ionic 
vibrational motion that is luq — uj c — u, in interaction picture and rotating wave approximation is given by jg, 



Hi = hfl[a + C>Q + ct-Oq ] + % r) c o+bd\a + h.c 



(1) 



Here a^(a) and b^(b) are creation(destruction) operators for vibrational phonon and cavity field photon, respectively, 
and v is trap frequency. The ion phonon and ion-cavity coupling constants are Q and <?, whereas ak(k = z, +, — ) are 
the Pauli operators qualifying the internal state of the ion. The operator Ok is defined as 



O k 



exp 



E 



^ P^P + kY- ' 



(2) 



with Lamb-Dicke (LD) parameters relative to the laser field and the cavity field denoted by 77 = r\i, and r\ = r) c 
respectively. In the limiting case t?l <S 1 and rj c <SC 1, Ofc=o,i — * 1, as such the relevant part of interaction picture 
Hamiltonian reduces to 



Hi = hfl[a + + er_] + hgrj c a+ba + <j-\r 



(3) 
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To obtain unitary time evolution of the system we work in the basis, \g, m, n), |e, m, n), \g, m — l,n — 1) , and 
|e, m — l,n — 1). Eigen values and eigen vectors of Hj are obtained analytically. Starting from a given initial state 
of the system, written in terms of eigen functions of Hi, the pure state ^(t) for the composite system is obtained by 
solving the time dependent Schrodinger equation, 

H I *(t)=ihj*(t). (4) 

The density operator for the composite system defined as p ABC (t) — \^f(t)) {9(t)\, can be used to obtain reduced 
density operators p A (t),p B (t), and fP (t) for subsystems A, B, and C, respectively. 



II. INITIAL STATES OF THE SYSTEM AND GHZ STATE GENERATION 

Consider that an ion is prepared initially in state 

*(0) = {cos{9) \g) + sin{9) |e)) |m = 0, n = 0) , (5) 

with the center of mass in lowest energy trap level (to = 0), cavity in it's vacuum state (n = 0), and < 9 < tt. 
Defining, a = \grj c , and /i = s/a 2 + Q 2 , we solve Eq.Q and verify that for interaction time t p such that \it p = pn, 
p=l,2,..., 

*(i p ) = [cos(6>) cos(aip) |g,0,0) + sin(<9) cos(ai p ) |e,0,0) 

—i sm(9) sin(atp) \g, 1,1) — i cos(0) sm(at p ) |e, 1, 1)] . (6) 

For the choice \i = 4a, 9 = qir(q = 0, 1) at instant t\ — [n/ (4a)] the system is found to be in state 

y( tl ) = ^-j=-(\g,Q,0)-i\e,l,l)), (7) 



which is a ma" : 




FIG. 1: Contour plots of negativity N B (= N° and linear entropy Sp (= Sf) as a function of # and scaled time T(= at), for 
subsystems B and C. Initial states are of type (i) that is separable at T = 0. Variables 9 and T are in degrees. 
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To get further insight into tripartite entanglement, the time evolution of entanglement for the following sets of 
initial states is investigated: 

(i) *(0) = {cos{9) \g) + sin{9) |e)) |0, 0) , < 9 < n, a separable state. 

(ii) *(0) = {cos{9) \g, 1) + sin{9) |e,0)) |0) , < 9 < w, cavity state is separable at t = 0. 

(iii) ^E'(O) = {cos{9) \g,/3) + sin{6) \e, —(3)) |0) , < 9 < n, cavity state is separable at t = 0, with the ion prepared 
in a Schrodinger cat state. The coherent state |/3) , is given by 

|#=exp(-tf) f;^L|m), (8) 

\ / m=0 * 

where \j3\ is the average number of vibrational quanta associated with center of mass motion. 

III. MEASURES OF ENTANGLEMENT AND PURITY OF A STATE 

For a given density matrix operator p ABC (t) at instant t, acting on composite Hilbert space associated with system 
{ABC), reduced density operator for subsystem A is defined as, p {t) — TrBc{p ABC '(*))■ The partial transpose of 
density operator p A with respect to subsystem A reads as 

cLa d B do 

P Ta = ^2 2222( i > m ' n \p A \ i; r ; s ) I 3, m , n ) ( i,r,s\ . (9) 
Negativity for subsystem A defined in terms of trace norm of partial transpose of density matrix as 

1 1 n T A 111 

N A = \\P_J (1Q) 

is equal to the modulus of sum of negative eigenvalues of operator p TA . N A is a measure of entanglement of subsystem 

4 with it's pnmnlpmpnf. ( FZCJ\ in thp rnrrmnsitp nnantiim «v«t.pm ( A T^C1\ A/Tpasnrps nf pnt.a.ncrlpTriPTit N an H N c For 




FIG. 2: Si, versus scaled time T{— at), for initial states of the type (i) separable at T = 0, (ii) cavity state separable at T = 0, 
and (iii) cavity state separable and ion in a Schrodinger cat state at T — 0, for the choice 9 — 90°, (3 = 1. 
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FIG. 3: Contour plots of linear entropy Sf 4 , as a function of 9 and T(= at) , for initial states of the types (i) separable at 
T — 0, (ii) cavity state separable at T — 0, and (iii) cavity state separable and ion in a Schrodinger cat state ((3 = 1) at T = 0. 

Although matrix p ABC (t) represents a pure state, the reduced density matrices, p A (t),p B (t) and p c (t) do not 
necessarily do so. Linear Entropy, Si defined as 

s=^4i(i-t7V)). (ii) 

is used as a measure of purity of a state. For a pure state Si — while for a maximally mixed state Si = 1. Here d is 
the number of modes available to a subsystem. 

IV. RESULTS AND CONCLUSIONS 

For the three sets of states described in section (2), we have obtained analytic expressions for ^ ABC (t). These 
expressions have been used for numerical calculations of entanglement measures N A ,N B and N c , for the system 
parameter choice of p,/a = 4. We have also calculated the reduced density matrices, p A (t),p B (t), p c (t) and from 
these the linear entropies S A , S B and Sf . Some of the results are displayed in Figs. (1) to (3). Fig. (1) shows contour 
plots of negativity N B (= N c ) and linear entropy S B (= Sf) as a function of 9(0 < 9 < n). and scaled time T(= at). 
Initial states are separable at T = (type (i)). A comparison of contour shapes denounces the fact that negativity 
and linear entropy do not produce the same ordering of states. However at zero and maximum entanglement, the 
states correspond to pure states and maximally mixed states, respectively. We find the subsystems A and B to be 
maximally entangled and maximally mixed at T = 45° and 135° for all possible choices of 9. 

Fig (2) displays for the choice 9 = 90°, linear entropy Si versus scaled time T(= at) in degrees, for subsystems A, 
B and C. Parts (i), (ii) and (iii) in the figure refer to initial states of the type (i), (ii), and (iii), as in section (2). 
For initial state (i), all three subsystems are found to be in maximally mixed state at T — 45°and 135°, hence are 
maximally entangled. For initial state (ii), although subsystems A and C become maximally entangled at T = 45°and 
135°, system B fails to reach maximum entanglement. For case (iii), when ion's vibrational state is a coherent state 
(/? = 1) at t = 0, entanglement is seen to increase with time for all the three subsystems. As expected, the inial state 
temperature has strong influence on generation of maximally entangled tripartite state. 

Contour plots of linear entropy Sf, as a function of and T(= at) , for initial states of the types (i), (ii), and (iii) 
are displayed in Fig. (3). For initial states (i) subsystem B is a qubit, for initial states (ii) subsystem B is a qutrit, 
and for initial states (iii) a large number of modes are available to system B, while A and C remain as two mode 



■5 



systems. A comparison of parts (i), (ii), and (iii) of Fig. (3) shows that for a given choice of 9, the time evolution 
of purity and entanglement of subsystem A depends strongly on the number of modes available to subsystem B. 
Time evolution of Si for system A in Figure (2) is a particular case of Fig. (3). We conclude that (a) the number 
of modes available to a subsystem not only determines the maximum entanglement of that subsystem but also has a 
strong influence on entanglement of other subsystems, b) at entanglement maxima and minima, linear entropy and 
negativity uniquely determine the nature of state, but the two measures do not induce the same ordering of states, 
c) for a special choice of system parameters maximally entangled tripartite two mode GHZ state is generated. The 
scheme presented for GHZ state generation is a single step process and is reduction free. 
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